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ABSTRACT 

This work presents a systematic investigation of the influence of weather conditions on the 
calibration errors by using Gaussian fitness, least chi-square linear fitness and wavelet trans¬ 
form to analyze the calibration coefficients from observations of the Ghinese Solar Broadband 
Radio Spectrometers (at frequency bands of 1.0-2.0 GHz, 2.6-3.8 GHz, and 5.2-7.6 GHz) during 
1997-2007. We found that calibration coefficients are influenced by the local air temperature. 
Considering the temperature correction, the calibration error will reduce by about 10% — 20% at 
2800 MHz. Based on the above investigation and the calibration corrections, we further study 
the radio emission of the quiet-Sun by using an appropriate hybrid model of the quiet-Sun at¬ 
mosphere. The results indicate that the numerical flux of the hybrid model is much closer to the 
observation flux than that of other ones. 

Subject headings: methods:data analysis - methods: numerical - Sumatmosphere - Sun: radio 
radiation 


1. Introduction 


Broadband radio spectrometers play an important role in observing and revealing the physical processes 
of the solar atmosphere, solar flares, coronal mass ejections, and other solar activities. When we utilize 
the data observed by spectrometers to study solar problems, calibration is a key procedure that dominates 
the reliability of the observational results. At present, there are several sola r broadband radio spectrome¬ 
ters running in the world, such as Phoenix at ETH Zur ich (100-4000 MHz ; Benz et al. ( 1991[B . Ondrejov 


Radiospectrograph i n the Czech Republic (800-5000 MHz; Ijiricka et al.l (|l993f) b Brazil Broadband Spectrom¬ 
eter (200-2500 MHz; Sawant et al. ( 200ll)i and th e Chin e se Solar Broadb and Radio Spectrometers (SBRS, 


1.10-2.06, 2.60-3.80, and 5.20-7.60 GHz; IFu et al.l (|l995[l . I.Ti et al.l (j2005[ B at the Huairou Solar Observing 


Station. All the above spectrometers are single-dish telescopes that receive the radio emission of the full 
solar disk without spatial resolutions. So far, calibration techniques of single-dish radio telescopes can be 
classified as follows: 


(1) Absolute calibration. Skillful practices are needed to determine the absolute calibration values by 
using an approximate gain of the antenna based on radar and communication techniques, together with an 
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approximate sensitivity derived from experiments with limit e d accura c y. It is still a difficu l t prob lem and 
usually only applied to polarimeters (IBroten fc Medd (119601) : iFindlavl (119661) : iTanaka et alJ (119731) : etc). 


(2) Relative calibration. This method needs reference emission sources to calculate the emission flux (or 
brightness temperature). It is usually applied to spectrometers since th e absolute calibration of spectrometers 
would be a huge and very complex mission at broadband frequencies ( Messmer et al. 19991) . 


(3) Nonlinear calibration. When the gain factor of the receiv er exceeds its no rmal range during strong 
solar radio burst, it is necessary to consider nonlinear calibration (|Yan et al.ll2002i) . 


Recently, we found that the cal ibration data ha s a strong correlation with the local air temperature at 
frequency of 2.6 - 3.8 GHz of SBRS (|Tan et al.ll2009l) . The local air temperature is an important factor that 
needs to be taken into account in the calibration procedure. We also analyze the impact of the Sun elevation 
and other weather conditions on the calibration. 


Moreover, the calibration of solar radio spectrometers is very important when we compare the observa¬ 
tions with theoretical work in the study of the quiet-Sun radio emission. It supplies an accurate flux spectrum, 
which help us to study the basic nature of the quiet-Sun emission and provi de a fundament al knowledge of 
the solar atmosphere. The radio quiet-Sun has been studied for about 70 yr. iMartvnI (|l946l) considered the 
quiet-Sun as a black body radiation and studied its radio spectrum. He showed the variation of radio bright- 
ness across t he solar disk at various frequencies and indicated that limb brightening should be observed. 
Smerd ( 1950l) studied the radio radiation fro m the quiet-Sun with nume rical analysis by applying radiative 


transfer equations to a typical ray trajectory ( Jaeger fc Westfold 1950l) . He presented a complete analysis 
of the quiet-Sun radio radiation, and the res ult is consistent with the observations at wavelengths of 3 cm , 


10 cm, 25 cm, 50 cm, 60 cm, 1.5 m & 3.5 m ( Pawsev fc Yabslev ( 19491) : 


Since then, more and more observatio ns and studies have been published (jAllenI (jl957^ ; I Tanaka et al.l (|l973l ) ; 


Christiansen fc Hindman (11951 )). 


Bastian et al. ( 19961 )). Kjmdu (1965) review ed the basic theory of radiative transform and propagation in 


solar radio astrophysics. I Tanaka et al.l (|l973l ) used the average value of the daily noon flux to obtain the radio 
flux spectrum of the quiet-Sun. Other works (INelson et al.l (|l985l) : IZirin et al.l (|l99ll) ) used the flux density 
of the quiet-Sun around the minimum of the solar cycle. The numerica l computation of the q ui et-Sun radio 
emission has been discussed in those papers (Bulk ( 1985 ): BenzI ( 1993 ): Selhorst et al. ( 2005 )). Benj ( 20091) 
pointed out that the radio emission of the quiet-Sun is a well-defined minimum radiation level when the Sun 
has no sunspots for s ome weeks. The presen ce of sunspots enhances the radio emission and produces a slowly 
variable component. IShibasaki et al.l (120111) reviewed radio emission of the quiet-Sun in observations and nu¬ 
merical computation. However, until now few works have paid attention to the study of the totally quiet-Sun 
radio flux spectrum at centimeter to meter wavelengths. Our work will contribute new observations and com¬ 
putations on this topic by improving the study in several aspects: (1) adopt new flux spectrum observations 
during solar cycle 23; (2) select the daily noon flux without sunspots ±3 days as the flux of the quiet-Sun; 
(3) utilize numerical integration inste ad of an approximate analytical function in n u merical computation o f 
radio emission; (4) use three models (IVernazza et al.l (|198ll) : ISelhorst et al.l (120051) : iFontenla et al.l (120111) ) 
of the solar atmosphere for comparison of numerical computations. We do not consider the magnetic field 
B and the scattering effect since the magne tic field B is weak i n the quiet-Sun and the scatterin g effect can 
be neglected at shorter meter wavelengths (|Aubier et al.l (|l97ll) : iTheiappa fc MacDowall I (120081) ). 


SBRS has observed the sun since 1994 and accumulated observation data for more than one solar 
cycle. In this work, we take the observation data of SBRS to analyze the daily calibration and increase 
the accuracy of the calibration procedure in section 2. Based on the radio flux spectrum of the quiet-Sun, 
section 3 presents numerical analysis of the quiet-Sun radio emission with improved numerical treatment 




























































































































- 3 - 


and several new models of the solar atmosphere. Section 4 gives the main conclusion and some discussions. 


2. Observations and Calibration 


Three spectrometers of SBRS (Fu et al. 1995, 2004; Ji et al. ( 2005ll l are located at the Huairou Solar 
Observing Station in China. The first spectrometer (SBRSl) was upgraded three times. SBRSl observed 
the Sun and only saved the burst data in the early years. It started routine and calibration observations in 
1999 October . It was under construction after 2006 July and started to work again in 2013. The second 
spectrometer (SBRS2) began to observe the Sun in 1996 September. The third spectrometer (SBRS3) began 
observing the Sun in 1999 August. SBRS2 and SBRS3 shared the same antenna, with a diameter of 3.2 
meter. The main information and performed parameters of SBRS are listed in table [T] During solar cycle 
23, the SBRS observed the Sun and did daily calibration at each noon (avoiding the effect of radio bursts), 
except for the instruments under maintenance. They supply plenty of observations and calibration data to 
analyze the quiet-Sun radio emission. In this work, we do not consider the daily calibration data in less than 
1 yr of observations because small data sets may result in larger error. We only select the daily calibration 
data from October 1999 to October 2004 of SBRSl, from 1997 to 2007 of SBRS2, and from 2000 to 2007 of 
SBRS3 except for some big data gaps (> lOdays) or without a noise source or termination (Left panel of 
FiglD). We go on analyze the daily calibration data with the daily noon flux from the National Geophysical 
Data Center (NGDC), and we study some treatments to reduce the calibration errors. At last, the flux 
spectrum of the quiet-Sun can be obtained from the observing data with calibration. 


Table 1: The main information and performance parameters of the component spectrometers 


Band 

(GHz) 

Dm 

(m) 

fre 

(MHz) 

Gadence 

(ms) 

Pol 

Ghan 

Observation Period 

SBRSl (1.0-2.0) 

7.3 

20.0 

100 

<10% 

50 

1994-2002 Jan 

SBRSl (1.10-2.06) 

7.3 

4.0 

5 

<10% 

240 

2002 May-2004 Oct 25 

SBRSl (1.10-1.34) 

7.3 

4.0 

1.25 

<10% 

240 

2004 Oct 26-2006 Jul 

SBRSl (1.10-2.10) 

7.4 

-2.78 

5 

<10% 

360 

2013 Aug 20-now 

SBRS2 (2.6-3.8) 

3.2 

10.0 

8 or 2* * 

<10% 

120 

1996 Sep-now 

SBRS3 (5.2-7.6) 

3.2 

20.0 

5 or 1.25* 

<10% 

120 

1999 Aug-now 


Note. — The first column ’Band’ is the frequency band of the spectrometer. ’Dm’ is the diameter of the antenna, ’/re’ is 
the frequency resolution. ’Cadence’ is the time cadence. ’Pol’ is the circular polarization accuracy. ’Chan’ is channel number. 

*the temporal resolution when the spectrometer works at a quarter of frequency band. 


2.1. Fundamental of Calibration 


Since our data are obtained by broadband spectrometers, this w ork mainly focu s on t he re lative cali¬ 
bratio n. The basic principle of the relative calibration was discussed in Messmer et al. ( 1999ll and Yan et al 


(|2002f) in detail. SBRS records a set of calibration data in daily routine observations. The left panel of FiglT] 
presents the schematic chart of the antenna and receiver system. The switch may turn to the antenna, noise 
source, and termination, respectively. Usually, the termination is a 50 ohm resistor. For each case there is 























- 4 - 
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Fig. 1.— Left panel is the schematic chart of the antenna and receiver system. Right panel is a figuration 
of the daily calibration data. 


an input to the receiver system. Ta is equivalent antenna temperature (Tq = Tgun or Ta = Tsky when the 
antenna points to the Sun or the sky, respectively). T„ and Ti are the equivalent temperature of the noise 
source and the terminal source, respectively. Therefore, there are four outputs at the computer at various 
frequencies: the Sun (marked as Rsun), sky background (Rsky), noise source (i?n), and terminal source 
They are showed in the right pan el of Fig U by the data processing software. When the gain factor (Gr) of 
the receiver is in the linear range (|Yan et al.ll2002r) . they comply with the following relationships: 


Rsun{v) = {Tsun + ■ Gr{v) 

(1) 

Rsky{v) = {Tsky +Tr{i^)) ■ Gr{iy) 

(2) 

Rn{v) = {Tn{v) + Tr{v)) ■ Gr{v) 

(3) 

Rt{v) = {Tt{v) + Tr{v)) ■ Gr{v) 

(4) 


Here v i s the observed frequency , an d Tr is equivalent temp erature of the receiver system, which is very 
small. From iBroten fc Medd (1960) and iMessmer et al.l ( 19991) . we have. 


Rsuniy') — 


‘2k B ■ Tsunjl^) 
Aeiy) 


(5) 


Fsun and Fsky are the radio fluxes of the quiet-Sun and sky, respectively. The effective area of the 
antenna aperture A^iy) is changeless and can be considered as a constant, ks is Boltzmann’s constant. 
The real flux of the sun Fq{v) should subtract the contribution of the sky Fsky{i^), i-e., equation ©• From 
equation m and © we have: 


Tb(^) — Fsuniy') Fskyiy') 


‘2k B Rsuniy') Rskyiy^} 

Af,{v) Gr{y) 


( 6 ) 


Gr{v) can be deduced with equation (|31) and (H]). Then equation ([5]) can be transformed into 
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Fo{y) 


2kB ■ {Tnjv) - Tt{v)) Rsun{l^) - Rsky{v) 
Aeiy) Rn{v) - Rt{v) 


(7) 


where Rsun{v), Rsky(y), Rn{v) and Rt(y) are daily calibration data. As written previously, ks and 
Ae(iy) are constant. T„ and Tt are equivalent temperatures of the noise source and the terminal source, 
respectively. The noise source and the terminal source are fixed electronic apparatus and usually stable 
under the steady environment and no interference. Thus, in equation ([8]), C(y) will be also stable under the 
steady environment and no interference. It is defined as calibration coefficient 


C{v) = 2kB ■ 


Ae{v) 


Fo{v) 


Rn{y) - Rtiy) 

Rsuni^^ Rskyiy') 


( 8 ) 


For any observed data Rsuniv) , the corresponding real flux Fq{v) of the Sun is calibrated as follows: 


Fq{^') — Fguni^iy^ Fgkyiy^ 


Rsuniy^ Rskyiy^ 
Rniv) - Rtiy) 


C{v) 


(9) 


Equation (jS]) is very convenient to do calibration for it does not require that we know the values of 
Tn{v) — Tt{v) and Af-iv). The bandpass flatness of the spectrum is usually not good for two reasons: 1) 
the gain factor Gr(i^) of the receiver varies along frequency; and 2) the antenna system Aeiv) will also vary 
more or less along frequency. The calibration with equation Q will eliminate the frequency property of the 
gain factor and flat the frequency property of the antenna system with Civ) of each frequency. The most 
important of calibration work is to analyze a suitable calibration coefficient Civ) with the right-hand side of 
expression ([5]) and reduce the calibration error. Usually we decide the daily noon flux from NGDC as Foiv). 
It can be downloaded from the NGDC Web sit^, which provides standard flux of solar radio emission at 
nine fixed frequencies (245, 410, 610, 1415, 2695, 2800, 4995, 8800, 15400MHz). We can calculate flux at 
any frequency between 245 and 15,400 MHz in using of linear interpolation. 


2.2. Impact of the weather in the calibration 

Generally, a constant coefficient is adopted in calibration. In practice, the calibration is influenced by 
the local air temperature more or less. Panels (b) and (c) of Figl5] show the comparison between the local 
air temperature and the calibration result at 2800 MHz. Panels (d) and (e) of Fig H] show that Rsunii'), 
Rskyiv)^ Rni^), and Rtiv) are all correlated with the local air temperature. From equation ([T]) and dD) we 
can deduce that the gain factor Gr(i/) is influenced by the local air temperature since Tgunit^) and Tgkyiv) 
have no relationship with air temperature and Triv) is very small. In fact, the gain factor Criv) of all bands 
of SBRS is influenced by the local air temperature more or less. Panels (d) and (e) also show two major 
jumps (arrows in the figure) at 2001 April 19 and 2002 November 09, respectively. We check the daybook 
of the observer and find that there is a change to the attenuation of the instrument. Panel (f) shows that 
the daily calibration coefficients Cd are correlated with the local air temperature. The bottom three panels 
of Figl2] indicate that the correlation between Cd and the local air temperature varied because adjustment 
of the instrument. We can deduce that Tniv) — Ttiv) in equation (|S]) is also influenced by the local air 


^http: //www.ngdc.noaa.gov/stp/SOLAR/ftpsolarradio.html 
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temperature. All these results imply that the electronic apparatuses of the instruments are influenced by the 
local air temperature. The relationship between the calibration result (panel (c) of Figj^]) and the humidity 
(panel (a) of Figl^]) or other weather conditions is not clear. The observation must be stopped when some 
extreme weather events occur (such as thunder, windstorms, rainstorms and heavy snow). Until now we 
have found no distinct evidence that the calibration was influenced by normal weather conditions except for 
air temperature. 









10 20 30 

Temperatrue °C 


0 10 20 
Temperatrue “C 


0 10 20 
Temperatrue “C 


Fig. 2.— (a and b) Local humidity and local air temperature, respectively; (c) difference between the 
daily noon flux of SBRS (calibrated with C) and the daily noon flux from NGDC; (d) Rsun (red small 
cross) and Rsky (blue dotted); (e) (red small cross) and Rt (blue dotted); (f) daily calibration coefficients 
Cd = Fq ■ ■ The bottom three panels plot the daily calibration coefficients vs. local air temperature 

for three terms. The sloping line is the least chi-square linear fitness of Cd and the local air temperature. 
All the radio data are at 2800 MHz. 


We will explain that the Sun elevation has a small effect on the calibration of this work. iTsuchiva fc MacDowall 


(|l965ll pointed out that solar flux can be measured without considering the weather condition (atmospheric 
absorption) at frequency lower than 17 GHz. Gonsidering the atmospheric absorption, the antenna temper¬ 
ature is given by equation as follow. 
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Ta — Tsource^ + Tsfey (10) 

T,ky= r T{T)e-^dT (11) 

Jo 

Tsource IS the real temperature at the emission source. T{t) is the temperature along the line of sight. 
T is the optical depth of the atmospheric absorption. 


T = 



K-dh-sec{z) 


( 12 ) 


where k is the absorption coefficient of the atmosphere. It is small (< 0.006 //cto) at frequencies lower 
than 17 GHz (|Tsuchiva fc MacDowall II1965I) . Moreover, the density of the water vapor and oxygen decreases 
rapidly along the height of atmosphere < IQkm above the ground level. The absorption k at the height of 
^ \Qkm will be about 1% of that of the ground level. Here h is height and z is the zenith distance. So 
the Sun elevation is equal to 90° — z. The daily calibration observation is done around the noon. The Sun 
elevation at noon is larger than 26°, considering that the latitude of Beijing is 40°. The sec(z) is smaller 
than 5.8 when the Sun elevation is greater than 10°. Thus the optical depth in equation m is also very 
small. To < 0.02. It has only a small effect on the temperature of the emission source. 


We also find no distinct evidence that the calibration was influenced by the Sun elevation. If the Sun 
elevation has a considerable effect on the observation, the following will result: (1) the variation profile of 
Rsun{oi Rsky) will show Considerable difference from that of i?„ (or Rt). The left panel of Fig[T]shows that 
Rn and Rt have no relationship with Sun elevation for no connection with antenna; (2) i?s„„(or Rsky) in 
summer (high elevation) will be larger than that in winter (low elevation). This is opposite to the observation 
shown in panel (d) of Figl2l From all the above analyses, the Sun elevation does not need to be considered 
in calibration. 


We need to understand two problems: (1) the variety of daily calibration data and (2) the relationship 
between the daily calibration coefficients and the local air temperature. The first is related to the daily 
calibration data and daily standard flux of solar radio emission. The daily calibration data was recorded 
by SBRS. They consist of 4 recorded parameters: Rsun, Rsky, Rn, and Rt (right panel of Fig[T|). The daily 
noon flux of radio emission at nine fixed frequencies can be downloaded from the NGDC Web site. Usually 
the observations from Learmonth Observatory are considered since its recording time (05:00 UTC) is very 
closed to that of SBRS (04:00 UTC). The second problem is related to the local air temperature data of 
Beijing, which are collected from China Meteorological Administration. 

The daily calibration coefficients are calculated by the right-hand side of equation dH) with the daily 
calibration data of SBRS during 1997 - 2007. There are several steps of pretreatment to exclude the abnormal 
data: (1) Look up the notebook of the observer, check out the day when the instrument was not in good 
operation (testing, maintenance, no noise source, etc.), and rule out the data during these cases; (2) get rid 
of the calibration data with interference, ±NaN value or minus value. The data are wrong and will result 
in unpredictable errors in the mathematical processing. The main reason for the ± NaN, minus, and big 
value of the daily calibration coefficient comes from the electromagnetism interference of the environment. 
So Rsky will be too big or even larger than Rsun- The relationship between the abnormal data and bad 
weather conditions is still not clear. We will list all the possible reasons in section 4. 
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Fig. 3.— Left: comparison between the average value C (long-dashed line) and the gaussian fitness value C 
(dashed line). Right: difference between the standard deviation of calibration of C and that of C. 


2.3. Analysis of the Calibration 

We take an example to analyze the calibration at 2800 MHz. At first, we compare the average value 
of the daily calibration coefficients (C) with the Gaussian fitness value of the daily calibration coefficients 
(C) at each frequency. They are C'(:/) in equation ([5]) when doing calibration. The Gaussian fitness value 
can exclude the contribution of big or very small values, which usually indicate radio bursts, sunspots or 
abnormal observations. The left panel of Figl3] shows the comparison between C and C at 2800 MHz. The 
vertical long-dashed line indicates values of C. The black solid line is the histogram of the daily calibration 
coefficient. The dashed line is the Gaussian fitness of the black solid line. The Gaussian fitness is: 

f{x) = AQe^ As-\-A 4 ^x-\-Abx‘^, z = ^——- (13) 

A2 

The X value of the maximum f{x) is C marked as vertical dashed line in FiglSl When \C-C\ < 0.02*17, 
the average value is close to the gaussian fitness value. This indicates that the observation of this frequency 
is reliable. In fact, C and C are close for most of the frequencies. The most probable value (MPV) is the 
maximum of the histogram (marked as star in FiglH]). The rms value is the square root of Cd (marked as 
cross in FiglSj. Both the calibration error of the MPV and the rms are lager than those of C and C. The 
right panel of Figl3] shows the comparison between the standard deviation {uav) of calibration with C and 
the standard deviation ( ctg ) of calibration with C . The standard deviation of calibration is calculated by 
equation m- 


a = stddev[FH — Fjv) (14) 

Here Fh is the daily noon flux calculated by equation © with the daily calibration data of SBRS and 
different C (C, C, etc.). Fjq is the daily noon flux from NGDG. For 82% of frequencies, (Jav are greater 
than ac- For the remaining 18% of frequencies, aav are less than or equal to (Jg because the gaussian fitness 
deviates from the histogram of the daily calibration coefficients. Thus, C is used as the constant coefficient 
of calibration. The calibration with a constant coefficient (C or C) is named the constant calibration. 
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Wavelet spectrum (Power) of daily calibration coefficient at 2800 MHz 
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Fig. 4.— Top: wavelet spectrum of daily calibration coefficients at 2800 MHz. The x-axis is the time (year), 
and y-axis is the period (days). The white cross is the daily calibration coefficients. The black smooth curve 
is the calibration coefficients after wavelet filtering. The black solid line shows the period of 365 days. The 
short white line near the black line of 365 days is the location of maximum power. The data gap represents 
no observation or during maintenance. Bottom: local air temperature. 


The other two sets of calibration coefficients are related to the local air temperature, including calibration 
coefficients with temperature correction (TC) and temperature-wavelet correction (TWC). Figl4] shows a 
strong correlation between the daily calibration coefficients (white cross in top panel) at 2800 MHz and local 
air temperature (black cross in bottom panel) during 1997 - 2007. Here we partition the data into three 
terms (1997 - 1999, 2000 - 2004, 2005 - 2007) according to the data gap during which the spectrometers are 
under maintenance. Those points that exceed the normal range during maintenance should be excluded. 
Then we decide in equation (1151) the linear relationship between the daily calibration coefficients Cd{t) and 
the local air temperature Tair{t) for each term. The variable t is the time count in days. Ci and C 2 can be 
obtained by least chi-square linear fitness between Cd{t) and Tair(t): 


Cd{t) — C\ + C 2 * Tairit) 


(15) 


Here Ci is the first term of the linear fitness. In the second term of the linear fitness, C 2 is the correct 
factor that is related to the local air temperature. For each day, the calibration coefficient Ctc = Ci+C2*Tair 
can be corrected by the local air temperature and replace C{v) in equation (|9]). This is the calibration with 
temperature correction TC. 


TWC is the calibration coefficient corrected by temperature-wavelet analysis. The top panel of FiglU 
plots the power spectrum ( red color) of Cd(t) a fter wavelet transform for three terms. The wavelet transform 
method was discussed by ISvch fc YanI (120021) in detail. Their work used a complex basis based on the 
Morlets wavelet which is well localized in both the scale and frequency plane. This makes it possible to 
rapidly reconstruct the signal even in the presence of a singularity value or in the absence of data. We first 
transform Cd(t) by wavelet to the wavelet data, which are complex numbers and give the power information 
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along the time and period. The period of the maximum power is about 365 days, and the half-power beam 
widths (HPBWs) are in the range of 260^530 days. We go on to set the value of the wavelet data beyond 
the HPBW as zero (that is, we only keep the wavelet data within the HPBWs), and we perform wavelet 
transform back to new coefficients C{t). This numerical process is wavelet filtering, which filters out the 
information we do are not concerned with. The new coefficients C(t) are plotted as a black smooth curve 
through Cd{t) (white cross). In equation (IT6l) . Tairit) is the local air temperature after wavelet filtering. So, 
Cl and C 2 can be obtained by least chi-square linear fitness between C{t) and Tairit): 


C{t) = Cl+C2*Ta^r{t) (16) 

Ctwc = C 1 +C 2 * Tair can replace C(v) in equation ([9]) when doing calibration. This is the calibration 
with TWC. 



0 100 200 300 100 200 300 100 200 300 100 200 300 100 200 300 


Fig. 5.— Each panel plots the daily noon flux of SBRS after C calibration vs. daily noon flux of NGDC. 
They are plotted as black crosses. The top three panels are the comparison of calibration results with C, 
Ctc and Ctwc at 2800 MHz. The bottom 11 panels are the yearly calibration results with C, Ctc and 
Ctwc at 2800 MHz. 

We compare the calibration results of three sets of calibration coefficients {C, Ctc and Ctwc) at 2800 
MHz in both the long term (more than 2 yr) and short term (1 yr). The standard deviations trc, ore, and 
ctwc are of calibration with C, Ctc and Ctwc, respectively. They are calculated by equation da. Each 
panel of FigIS] plots the daily noon flux observed by SBRS (calibration with C) versus the flux of NGDC. 
The upper three panels are the calibration results for three long observation terms. The bottom 11 panels 
show the yearly calibration results. The cttc or ctwc are about 10% ^ 20% smaller than ctg at various 
observation terms. The difference between gtc and otwc is small, ^ 2%. In short observation terms, ctc 
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ls the smallest. While in long observation terms, cftwc is a little smaller than arc in most cases. We 
conclude that TC calibration is better for short observation terms, while TWC calibration is better for long 
observation terms. The relative standard deviations (RSDs) of calibration with Ctc and Ctwc are less 
than 10% in all cases, while RSDs of calibration with C are greater than 10% sometimes because it cannot 
correct the influence of the air temperature. The RSD of the calibration will be larger when the influence of 
air temperature is stronger. The Ctc of 2007 is used in the calibration during the years of 2008 and 2009. 
At 2800MHz, arc = 6.9s./.u is about the same as for the years of 2004 - 2007. In practice, the calibration 
coefficient C'(^) should be updated annually. 

The daily calibration data of SBRSl and SBRS3 also have a strong relationship with the local air 
temperature. The relationship between Cd and the local air temperature is quite complex because SBRSl 
and SBRS3 were examined and repaired many times. In short observing terms less than 1 yr, there is no 
significant difference (< 5%) between constant calibration and TC (or TWC) calibration. The constant 
coefficients are used in calibration of them. At the 1.0-2.0 GHz band, the daily noon flux of SBRSl at 1420 
MHz is compared with the daily noon flux of NGDC data at 1415 MHz. The standard deviations are varied 
within 9 ^ 12s ./.m in different observation terms. The same is done at 1.10-2.06 GHz band. The standard 
deviations of calibration are varied within 8 ^ 12s./.u at 1420 MHz in different observation terms. At the 
5.2-7.6 GHz band, there are no observation data from NGDG. The linear interpolated values between 4995 
and 8800 MHz of NGDG are used as comparison. The standard deviations of calibration are varied within 
5 ^ 18s./.u at 5900 MHz in different observation terms. The observation term with the smallest standard 
deviation indicates a stable system and few occasional interferences, therefore signifying the best status of 
the instrument. 


2.4. Observation Results 


The daily noon flux observed by SBRS is calibrated by equation ([9]) with the calibration coefficient C{i'). 
At the 2.6 - 3.8 GHz band, C{i') is Ctc- At other frequency bands, C{v) is the constant coefficient. The left 
panel of FiglUplots the daily noon flux at 1415, 2800 and 5900 MHz observed by SBRS during solar cycle 23. 
In each band of SBRS, we plot the same (or nearby) frequency as nine fixed frequencies of NGDC. The top 
right panel plots the correlation coefficients between the sunspot numbers and the daily noon flux at variou s 
frequencies. The correlation coefficients are consistent with previous works ( Christiansen fc Hindman Ifiiilli . 
The bottom right panel plots the sunspot numbers and the daily noon flux at 2800 MHz in 2008 August, 
as an example. There are still small sunspots or weak active regions during August 21-22 even in the solar 
minimum. In order to obtain the pure radio flux of the quiet-Sun, we exclude the daily noon flux with 
sunspots ±3 days (gray cross in bottom right panel of FiglHl). During the solar minimum period of 2006- 
2009, 317 days are selected as the candidates. The Gaussian fitness of the daily noon flux of these 317 days 
is decided as radio flux of the quiet-Sun statistically. They are 55, 67, and 68 sfu at 1415, 2695, and 2800 
MHz of SBRS, respectively. The corresponding radio flux of the quiet-Sun from NGDG are 11, 27, 35, 55, 
67, 69, 118, 220, and 519 sfu at the nine fixed frequencies, respectively. 


3. Theoretical Analysis of the Quiet-Sun Radio Emission 

The quiet-Sun radio emission is thermal radiation originated from the ambient plasma in absence of 
solar activity. The mechanism is bremsstrahlung from electrons interacting with ions in the presence of a 
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Daily noon flux at three frequencies of SBRS/Huairou 
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Fig. 6.— Left: daily noon flux at 1415, 2800, and 5900 MHz observed by SBRS in solar cycle 23. The data 
gap is without daily calibration or during maintenance. Top right panel: correlation coefficients between the 
sunspot numbers and the daily noon flux at various frequencies of NGDC and SBRS. Bottom right panel: 
sunspots numbers (star) and the daily noon flux (cross) at 2800 MHz in 2008 August as an example. The 
gray cross is excluded because there are sunspots ±3 days. 


relatively weak magnetic field (IShibasaki et al.N2011r) . In an irregular propagation medium, a wave cannot 
be represented by a single ray. Small fluctuations in density or magneti c field will distort the incident 
plane wave as the wave phase propagates at different speeds. I Bend (|l993l l pointed out that the evidence 
of scattering of solar radio emission is ambiguous, while the scattering hypothesis has bee n successfully 
verifie d for pulsars and irregula r ities o f the interstellar medium. Moreover, other papers (lAubier et al 


(1971); Theiaopa fc MacDowall ( 2008ll l concluded that (1) the scattering effect decreases the intensity of 
the radio emission and enlarges the size of the radio source, and (2) the scattering effect increases with 
the radio wavelength (oc A^) and can be neglected at shorter meter wavelengths. We mainly studied the 
quiet-Sun radio emission at the frequency range of 245 - 15400 MHz in this work; therefore, the scattering 
effect can be ignored. In nume rical analysis o f the quiet-S u n radio emiss ion, most previous works are based 
on the theoretical treatments ( Smerd ( 1950ll : Dulk ( 1985 1: Benz ( 1993ll l. The quiet-Sun radio emission is 
calculated by using equations and treatments as follows. 


The radiation transfer equation is 


_ V I 

ds ^2 


(17) 


Here I is the specific intensity of the radiation at frequency of zz; t] is the volume emissivity, k is the 
absorption coefficient, and /i is the refractive index of the medium. The optical depth and the path element 
have the relationship dr = nds. In this work, the electron temperature T can be treated as uniform in 
small segments because of the entirely numerical integration. Thus, under conditions of thermodynamic 
equilibrium in small segment, we have rj = fpKB{T). In the radio frequency band, the Rayleigh-Jeans 
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approximation is B{T) = Hence, the solution of equation (flTl) can be obtained: 


Here the intensity Jq and refractive index fj,Q are at an optical depth of tq. kb is the Boltzmann’s 
constant. 



p in units of Rq 


Fig. 7.— Ray trajectory of 245 MHz in polar coordinates. The bending curves of various colors indicate the 
ray trajectory of various distances b from the center Sun-Earth line. The normal line of the refraction at 
the turning point is plotted as a dashed line. The arrow indicate the propagation direction of the ray to the 
observer. 


The ray treatment of radiation ( Jaeger fc Westfold 1950[) is based on the refraction of the ray path in 
the solar atmosphere. The equation of the path and the absorption of the rays can be deduced from Snell’s 
law. FiglT] shows the ray trajectory in polar coordinates {p,0). Here p = R/Rq is the distance from one 
point of the ray path to the solar center. These rays will emerge from the solar atmosphere parallel to the 
observer at a distance b from the center Sun-Earth line. Both p and b are in units of the sun’s optical radius 
Rq (6.95 X lO^^cm). The refractive index p of an ionized medium decreases with increasing electron density. 
It follows that a ray passing through the solar atmosphere experiences continuous bending by refraction. 
The point where the direction of propagation changes from that of decreasing p to that of increasing p is 
referred to as the ’reflection point’ or, b etter, as the turning point. All these rays are calculated with the 
equation m (| Jaeger fc Westfolcnil950ll and the solar atmosphere model in Fontenla et al. (1993, 2011) 
paper. The path element ds of a trajectory b is given by equation (uni) 


Jp p^/p^p"^ - 62 


(19) 


RQdp 


ds = Rq \/ {dpY + {pdOy 


( 20 ) 
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The optical depth can be deduced from the equation of the path (l20l) and the absorption of the rays. 
The optical depth between two points (pi and / 02 ) is as follows: 


fP2 

TiAb) = / 

j Ol 


kR(7) 


PI V^l - 62/^V- 


:dp 


( 21 ) 


The refractive index is ~ 1 — • (1 ± ^|cos0|)“^ ~ 1 — because the ray frequency is much 

greater than the local gyrofrequency w ^ fie along its propagation in the solar atmosphere o f the qu i et-Su n 
(Benz 1993 !). The formulae o f absorption coefficient k d i ffer in different papers. Some papers ( Smerd ( 1950 1: 


Bracewell fc Preston ( 1956ll : Theiappa fc MacDowall ( 201 oIB used k deduced by classical collision theory. 


The classical collision theory defined the absorpt ion coefficien t k = Here i'e,i is the collision 


frequency of thermal electrons/ions. Some papers (IDulkl (119851): iGarv et al.l (|l990r) l used k with the Gaunt 
factor ded uced by quantum the o ry. T h e formula of n in Dj^ (119^) is the same as the approximate analytic 
formula of Novikov fc Thornel ( 1973[l: Rybicki fc Lightman (19^ ) at radi o wavelength {hv < k^T). We 
compare the Gaunt factor value of loulk (1985) with that of vanHoof et al. ( 2014 1 at the solar atmosphere 
from 100 Mhz to 30 GHz and find that the difference is of < 2% which has a very s mall impact on the 
numerical result. Thus in this work we still use the value calculated by the equation of iDulkl (jl985|). Table 
[2] lists the equations and parameters from different papers. They are all approximatively in the range of 
under the condition of the solar atmosphere since the classical collision theory is valid 


0.15 - 0.25- 




and close to the quantum theory when hv ^ ki^T. 


Table 2: Parameters of collision frequency Ve^i^ average Gaunt Factor g^f(y,T), and absorption coefficient 

K. 


1 

2 

3 

4 

Paper 

I'ep (Cgs) 


K (cgs) 

-—!-. rn — 


Bracewell (1956) 
Bulk (1985) 
Bulk (1985) 
Thejappa (2010) 




2 2 ^ 3/2 


^[18.2+ M^)] 

^[24.5+ Zn(^)] 


[18-2 + T < 2 X 10 

[24.5 + ln{^)] T > 2 X 10 ^ 

+ ^^(^T~)] (corona) 


2jql/3l J 
0.2N^ 

5 


4.36iV 

T3/2 


[10.81+ Zn(^)] 


From equation (I21|) and k in Table [21 we have 


Tl,2{b) = 


9.78 X lO-^Rrr, TT 


_ >2 g^^{v,T)N^ 

VsJp, 


dp 


( 22 ) 


N and T are the electron density and electron temperature in the solar atmosphere, respectively. They 
will be discussed in the next subsection. Equation (I22|) has a singular i ty point when u = b/p , p = p n- This 


point was named the turning point or ’reflection point’ ( SmerdI (1950); Bracewell fc Preston ( 19561 1) where 


the direction of propagation changes from that of decreasing p to that of increasing p. The appendix proves 
that the integration is convergent near the singularity point so long as the electron density N increased 
limitedly and monotonously along the height p. Therefore, when p > pQ, the integration can be calculated 
with numerical integration. 
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3.1. The Electron Density and Temperature 



Height above photosphere {km) 


Fig. 8.— The electron density (top panel) and electron temperat ure (b o ttom panel) distribution of different 
models along the height above photosphere. The classical model lAllenI (|l947T ) is plotted as long dashed line 
only at height of higher than lO'^km. 


There are many models of the electron density and electron temperature of the quiet-Sun atmosphere. 
The VAL III model (IVernazza et al.l 119811) determined semi-empirical models for six components of the 
quiet solar chromosphere in using of EUV observations. The series of FAL models (Fontenla et al. 1993, 
2009, and 2011) built the semi-empirical models with the optical continuum and EUV/FUV observation. 
These two models are excellent for reconstructing the optical and ultraviolet observations but still have 
some discre pancies when describ ing the radio observations. The discrepancies will be i llustrated in the next 
subsection. ISelhorst et al.l (|2005ll proposed a hybrid model that uses the FAL C model (IFontenla et al.lll993ll 
from th e phot osphere to 1800A:m, IZirin et al.l (119911) model from 1800 to 3500km in the chromosphere, and 
Gabrieli (1992) model from 3500 to dOOOOfcm in the transition region and corona. These three models are 


selected in this work because the numerical results are not far apart from the observations. Figl8] shows 
the electron densit y and electron t emper ature distribution of different models along the height above the 
photosphere. The IFontenla et al.l (120111) model is plotted as dashed line. It gave parameters from the 
chromosphere to the high c orona until the height of ^ 2 x lO^fcm. In the higher corona, the electron density 
is given by the lAllenI (|l947l) model as follows: 


N = 10®(1.55p-® -b 2.99/9 -^®).cto"3 (23) 

The Allen model (long-dashed line in FiglSj) will match the value of the Fontenla model at the height of 
^ 2 X lO^km. The Vernazza model is plotted as a cross in FiglS] It gave parameters in the chromosphere and 
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transit region until the height of 2439fcm. We decide to use the same parameters as the Fontenla model in the 
corona and the Allen model in the higher corona as before. They are plotted as gray crosses in the figure. The 
Selhorst model is plotted as a solid line in FiglS] It gives parameters from the chromosphere to the corona 
until the height of 3.93 x lO^fcm. We decide to use the Allen model after modification (gray long-dashed 
line in FiglS]) in the higher corona. Then three hybrid models that we used in this work are the FAL-I-Allen 
(F-l-A) model, VAL-I-FAL-I-Allen (V-l-F-l-A) model, and Selhorst -I- Allen (S-l-A) model. For all the models, 
the electron temperature in the higher corona is equal to the last value given by the corresponding model. 
We will compare the numerical results of three hybrid models with the observations in the subsection 3.3. 


3.2. The Optical Depth and Brightness Temperature across the Solar Disk 
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Fig. 9.— Left panel plots the absorption (long dashed line) and the optical depth r„_oo (real line) along the 
solar height, for each of the central line at eight fixed frequencies. Right panel is the brightne ss temperature 
across t he solar disk at eight fixed frequencies. Black lines are the numerical result with ISelhorst et al 


( 2005 )4- Allen ( 1947 ) model. Dashed lines are the results of approximation solution in SmerdT i 19501) . 


We do numerical integration entirely in this paper. In equation the optical depth r„_„+i(&) is 

calculated for about 500 points from the turning point (inner limit) to the point after which the contribution 
is very small (outer limit). The appendix proves that the integration of equation (1221) near the turning 
point is convergent. Beyond the outer limit, the absorption is very small (< 10“^° km~^ in this work) 
thus, the optical depth can be ignored. The 500 points are decided as follows: (1) the zero point is the 
turning point po; (2) assume the first point pi = po + Api, Api < li?©, and calculate the optical depth 
Top between the two points; (3) set the midpoint between point 0 and point 1, pm = Po + (pi — Po)/2, and 
calculate the optical depth of tq,™ and Tm,i, respectively; (4) if | '^o.m+Tm,i-To,i | ^ change the midpoint 

to point 1 and set a new midpoint. Do this circulatory calculation until | '^o,m+Tm,i-To^i | ^ iq-2 ^ 
calculation error is very small. Then the first point pi = p^ is decided. Usually, the first point is decided as 
Api = 10“^^ ^ 10“^°1?Q by experience. (5) Decide the rest points with this method. It should be taken care 
that the electron density and electron temperature vary abruptly in the transition region. The step should 
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be < Ap„ < 1.5 X Ap„_i to avoid the big error during calculation. For the points beyond the 

outer limit, the optical depths are approximated as zero. The left panel of FiglH] plots the absorption and 
optical depth t„^oo of the central line calculated with the Selhorst + Allen model at eight fixed frequencies 
as an example. The right panel of FiglH] shows the brightness temperature (or radiation intensity), which is 
calculated with equation (flSl) and (|2^ for different b at eight fixed frequencies. Blac k lines are the numerical 
result of the Selhorst + Allen model in this work. Dashed lines are the result of SmerdI ( 1950ll with the 
approximation solution. It shows differences between two results. 

Figini illustrates the main generating region of the quiet-Sun radio emissi on at v a ried f r equencies. The 


result is consistent in the main but a little different from the classical results ( Smerd ( 1950ll : DelaLuz et al 
( 20inl ll. At high frequency (3 ~ 30GHz), the absorption (< and optical depth (< 0.01) is 

very small in the corona. The emissivity rj = ii^kB{T) are very small in the corona. Thus, the ther¬ 
mal bremsstrahlung emission of the quiet-Sun at high frequency (3 — 30GHz) is mainly produced in the 
chromosphere. At low frequency (< 300MHz), the radio emission cannot propagate into the chromo¬ 
sphere. The thermal bremsstrahlung emission of the quiet-Sun at low frequency (< 300MHz) is mainly 
produced in the corona. The thermal bremsstrahlung emission of the quiet-Sun at intermediate frequency 
{300MHz < V < 3GHz) is produced in both the chromosphere and the corona. The limb brightening 
appears obviously at a frequency range of 200MHz < v < 30GHz. The brightness temperature of the solar 
center radio emission is consistent with the electron temperature of the main generating region. 


3.3. Brightness Temperature and Flux Density of the Quiet-Sun 




Fig. 10.— Left: brightness temp erature spectrum of the center line of the quiet-Sun. The errors bar are 


derived from the observations in Zirinetal 


(jl99ll l. Middle: comparison between numerical results and 
observations of flux spectrum of the quiet-Sun. Th e error bars are derived from observations of SBRS 
and NGDC, respectively. The observation result in Benz ( 2009ll has an accuracy of about 10%. Right: 
enlargement of a small rectangle from the middle panel. It shows the observation fluxs (thick blue lines) and 
error bars (thin blue lines) of TC and no TC in the range of 2500-4000 MHz as an example. 


























































- 18 - 


We calculate the brightness temperature spectrum and flux spectrum of the quiet-Sun radio emission 
for three hybrid models (F+A, V+F+A, S+A) at various frequencies. The total amount of radiation per 
unit time, unit frequency interval, and unit angle from the Sun to a distant observer is given by 


E = 2ttRI / I{d)d6d, 


(24) 


in units of erg.cm ‘‘Hz ^. In practice, the upper limit of this integral is replaced by a finite value dp, 


out of which the radiation contribution is very small. The flux density of the quiet-Sun can be transformed 
with the conversion Is.f.u = Hz~^ = 10~^^erg.cm~^Hz~^s~^. Then we compare the numerical 

results with observations and find the difference. The left panel of FigHUlshows the comparison of brightness 
temperature sp ectrum betw een the observat i ons a nd numerical results of different models. We chosse the 
observation s of iBena (|2009l) and IZirin et al.l (119911) as the standard spectrum because the flux spectrum of 
Benz ( 2 OO 9 I 1 is well consistent with the observations of SBRS-/F1SO and NGDC (middle panel of Fig 1^. A t 
frequency o f IbAO O MHz or nea rby, the numerical results of all models and all the obseryations (IFuerstl (Il980li ; 
Zirin et al. ( 1991 ): Benz ( 20091) 1 are consistent. This indicates that the parameters in the low chromosphere 
of all the models fit the obseryation well. At frequency of 3 — lOGHz, the numerical results of the F-(-A 
m odel (mainly F ontenla model) and V-|-F-(-A model (mainly Vernazza model) are close to the obseryations 


of 


Zirin et al 


Fuerst ( 1980 h while the numerical results of the Selhorst-I-Allen model are close to the obseryations of 


(1991) and Benz ( 20091) . The radio emission at the frequency of 300MHz — 3GHz comes from 


the transition region and corona. The numerical results of the Selhorst-|-Allen model (red solid line) are 
close to the observations of Benz (cyan solid line) well. Thus, we think that from the transition region to 
the low corona the parameters of the Selhorst model fit the observations well. It is complex to estimate 
wether the parameters in the higher corona are good or not owing to a lack of observations at low frequency 
(< 300MHz). 


Table 3: The center-line brightness temperature and flux density of the quiet-Sun at various frequencies 


v[MHz] 

Tc[K] 

F + A 

Tc[K] 

V + F + A 

Tc[K] 

S + A 

Tc[K] 

Benz 

Fq 

F + A 

Fq 

V + F + A 

Fe 

S + A 

Fq 

Benz 

Fq 

Fn 

Fq 

Fh 

245 

1.00 X 10^ 

1.00 X 10® 

9.69 X 10® 

7.88 X 10® 

17.3 

17.3 

15.7 

11.2 

11 


410 

7.93 X 10^ 

7.89 X 10® 

5.90 X 10® 

5.50 X 10® 

41.4 

21.3 

30.3 

22.2 

27 


610 

5.56 X 10^ 

5.27 X 10® 

2.84 X 10® 

3.58 X 10® 

69.8 

68.9 

41.9 

32.3 

35 


1415 

1.68 X 10® 

1.48 X 10® 

7.43 X 10* 

9.70 X 10* 

139 

127 

60.8 

46.9 

55 

55 

2695 

6.21 X 10'* 

5.85 X 10* 

2.93 X 10* 

4.19 X 10* 

185 

167 

77.9 

64.7 

67 

67 

2800 

5.76 X 10* 

5.49 X 10* 

2.76 X 10* 

3.68 X 10* 

189 

172 

80.3 

66.8 

69 

68 

4995 

2.74 X 10"* 

3.25 X 10* 

1.62 X 10* 

1.88 X 10* 

246 

251 

121 

107 

118 


8800 

1.49 X 10* 

2.10 X 10* 

1.22 X 10* 

1.38 X 10* 

359 

450 

233 

235 

220 


15400 

9.46 X 10® 

1.29 X 10* 

1.04 X 10* 

1.16 X 10* 

612 

832 

562 

599 

519 



Note. — The first column is frequency. The next three columns are numerical results of center-line brightness temperature 
w ith dif f erent models (FAL-hAllen; VAL+FAL+Allen; Selhorst+Allen). The fifth column is center-line brightness temperature 
in iBenj |200^. From sixth to eighth columns are nume rical results of flux density with different models. The last three column 
are observation results of flux density of the quiet-Sun teenj {200^; NGDC online data; and SBRS). 


The middle panel of Fig jlQI olots the flux spectrum of the quiet-Sun radio emission. The quiet-Sun fluxes 
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of SBRS observations (Gaussian fitness value) and the error bars (Scr) are plotted as blue crosses and short 
vertical lines. The quiet-Sun fluxes of NGDC (Gaussian fitness value) are plotted as black stars. The NGDG 
did not present the error of solar flux. Here we use three times the standard deviation (3 ct) of the quiet-Sun 
data from NGDC as error bars (black short vertical lines). The standard deviation is calculated separately 
at the low value part and high value part. We find that both the quiet-Sun flux spectra of SBRS and NGDC 
fit the observation of Benz (cy an sol i d lines ) well. Thus, our work agrees with the brightness temperature 
spectrum and flux spectrum of Benj ( 2009l l as the standard spectrum of the quiet-Sun radio emission. All 
the values of brightness temperature and flux density of the quiet-Sun at various frequencies are listed in 
table [21 The right panel of Figlin] shows the observation flux and error bars of TC (solid blue line) and 
no TC (dashed blue line) in the range of 2500 - 4000 MHz as an example. It indicates that calibration of 
TC is more accurate than that of no TC. The red lines are the same numerical results of solar models as 
in the middle panel of Fig {TUI The cyan line is the observation of Benz. As analyzed in subsection 2.3, the 
difference between TC and TWC calibration is very small (^ 2%). The observation flux and error bars of 
TWC are not plotted together. 


4. Conclusion and Discussion 


SBRS has been observing the Sun and obtaining plentiful data since 1994. This work adopted the 
observations of SBRS to investigate the calibration procedure and study the quiet-Sun radio emission. The 
calibration gives accurate observations which is basic and important in the study of the quiet-Sun radio 
emission, while the study of the quiet-Sun radio emission is a scientific extension of the former part. We 
first study various impacts to the calibration results and conclude that: 

1) Generally, the calibration coefficient is constant and should be upgraded annually. Actually, the 
calibration result with constant coefficient is found to be related with the local air temperature at all 
frequency bands of SBRS. One possible reason is that the electronic apparatus of the instruments are 
influenced by the local air temperature. Thus, the correlation between calibration and local air temperature 
will vary if there is an adjustment to the instrument. 

2) The relationship between calibration results (panel (c) of Fig|21) and the humidity (panel (a) of Fig|2]) 
or other weather conditions is not clear. There is no distinct evidence that the calibration was influenced by 
normal weather conditions except for air temperature. 


3) The Sun elevation has a small effect on the calibration of this w ork. The absorption of the atmo sphere 
should be considered only when the frequency is higher than 17 GHz ( Tsuchiva fc MacDowall 1965[) or the 
Sun elevation is lower than 10°. 


The calibration accuracy is also influenced by the occasional abnormal signal. Some possible reasons 
for the abnormal signal are listed as follows: 

1) Various kinds of interference: wireless communication, plane and airport, satellite, vehicle, lightning, 
etc. These will result in large values. 


2) Unstability of the feed, cable, or noise source. 


3) Sometimes, the gain factor may be out of normal range when there is a strong signal (|Yan et al 


20021 1. The beam of the antenna will offset the sun center if the tracking control is out of normal range. 


These will result in big or saturated values, or small values. The relationship between the abnormal signal 
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and bad weather conditions is still not clear. The calibration with equation @ will eliminate the frequency 
property of gain factor Gr{v), and flatten the frequency property of the antenna system with C{v) of each 
frequency. The calibration with equation Q does not need the values of Ae(y) and Tn{v) — Tt{v). 


In order to improve the calibration accuracy, we investigate the influence and properties of the instrument 
from the data analysis and comparison between different calibration coefficients. All the investigations are 
under the fundamental of calibration. First, we exclude the abnormal data that are not well observed or are 
influenced by the interference. Then we compare the calibration results of four sets of calibration coefficients, 
including average value (C), the Gaussian fitness value (C), constant coefficients after temperature correction 
(Ctc), and constant coefficients after temperature-wavelet correction (Ctwc)- The main analysis results are 
as follows: (1) In the 2.6 - 3.8 GHz band, the calibration errors ac are smaller than aav at 82% of frequencies. 
Comparing with other constant coefficients of calibration, C is the best. (2) At 2800 MHz, arc or utwc are 
about 10% — 20% smaller than ac at various observation terms. The RSDs of arc and arwc are less than 
10%, while RSDs of ac are greater than 10% sometimes because of the influence of the temperature. The 
Ctc is used in the calibration in the years of 2008 and 2009. The calibration error axe = 6-9s.f.u is about 
the same as that in the years of 2004 - 2007. (3) The calibrations of SBRSl and SBRS3 also have a strong 
relationship with the local air temperature. But there is no significant difference between TC (or TWC) 
calibration and constant calibration of SBRSl and SBRS3, which were examined and repaired many times. 
C is used in the calibration for several short observation terms. The standard deviations of calibration are 
varied within 8 — 12s./.u at 1415 MHz and within 5 — 18s./.u at 5900 MHz for different observation terms, 
respectively. The observation term with the smallest standard deviation indicates a stable system and few 
occasional interferences, therefore signifying the best status of the instrument. 


The daily noon fluxes of SBRS are calibrated with improved C{iy). We select the daily noon flux without 
sunspots ±3 days as the flux of the quiet-Sun during the solar minimum period of 2006 - 2009. The spectrum 
of daily noon flux observed by SBRS approaches the spectrum from NGDC we ll. The daily noon flux is 


corre lated with sunspot numbers. The result is consistent with previous works (jGhristiansen fc Hindman 


1951 11 ■ 


Based on the above investigations, we further study the radio emission from the quiet-Sun. The flux 
spectrum is strictly the observation without sunspots during the solar minimum period of 2006 - 2009. The 
numerical simulation of the quiet-Sun radio emission in this work is improved with several semi-empirical solar 
models and entirely numerical integration instead of approximate solution. The t heoretical and treatment ha s 
been testified well b e cause i t can deduce the same result as that of many papers (| Jaeger fc Westfold (Il950ll : 
Theiaopa fc Kundu ( 1992 1: Theiappa fc MacDowall ( 2010ll : etc). Different works have different treatments, 
such as (1) different approximations of radiation transformations, (2) different absorption coefficients k, (3) 
consideration of refractive index, scattering effect, or magnetic field. This work uses the radiation transfer 
equation and bremsstrahlung mechanism of emission without approximation and without scattering effect 
and magneti c field. The ray treatmen t of radiation is based on the refraction of the ray path in the solar 
atmosphere ([Jaeger fc Westfold 11195011 . The refractiv e index and absorption coefficient k injDul^ylWS^ 


are co nsidered since the Gaunt factor in iDulkl (jl985[ l is very close to the numerical value of IvanHoof et al 
(|2ni4ll . 


The VAL (jVernazza et al.lll98lll and FAL (jFontenla et al.ll201lll models are excellent for reconstructing 
the optical and ultr aviolet observations but still have some discrepancies in describing the radio observation. 
Zhang et al.l (j2n0l[ l discuss the discrepancy and attribute t he difference to an un derestimation of Fe abun¬ 
dance used in the calculation of the UV line emission. Thus, ISelhorst et al.l (j2005l l proposed a hybrid model, 
of which the transition region is about 1000 km higher than in typical VAL and FAL models. We regulate 
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solar models with a combination of different models and compare the numerical results with the observa¬ 
tions. The numerical results of the models should be within the confines (error bars) of observations. Finally, 
we find an appropriate hybrid model that uses the Selhorst model from the chromosphere to the height of 
3.93 X above the photosphere and the Allen model in the higher corona. The comparison (Fig fTUl) 

between the numerical results and the observations indicates the followings: (1) The numerical results of all 
models and all observations are consistent at 15400 MHz or nearby. The solar models in the low chromo¬ 
sphere are unanimous. The optical depth at 15400 MHz is very small in the transition region or higher no 
matter what model is selected. (2) The argument happened at the transition region or nearby, which is the 
main emission region at 300MHz — lOGHz. Only the flux spectrum of the Selhorst model is within the error 
bars of the observation. The theoretical results of the VAL and FAL models are about two times larger than 
observations. At the frequency of 300MHz — lOGHz the optical depth varied mainly within 0.01 < r < 10 
in the transition region or higher. It is difficult to modulate the solar models only with radio observations 
because each modulation in this region will impact on the results of the whole band of 300MHz — lOGHz. 
But the radio observations will help us to validate the models we known. (3) In the corona, the numerical 
results of the Allen model are close to the observations at the frequency of 100 — 300MHz. There is small 
argument on the solar models in the corona. But in this work we can not assure exactly the parameters in 
the corona as a result of few observations and the scattering effect in low frequency (< 300MHz). (4) The 
observation with higher accurate calibration will help us further qualify the empirical/semi-empirical solar 
models. The solar model can be decided when the theoretical result is within the error. 


The brightness temperature distribution across the solar disk in this work (right panel of Figl9 ) shows 
that the limb brightening happens in the frequency range of > 200MHz — 30GHz. iKundu et al.l (jl977n 
showed that no limb brightening happens at the frequency of lower than ^ 121MHz. Mercier & Chambe 
(2009, 2012) studied the radio images of the quiet-Sun at 150—450M7Jz but did not discuss limb brightening. 
However, their papers give us some clue that the limb brightening increased with the frequ ency in the range o f 
300—450MiJz. Some works also fou nd the limb brighteni ng at a high frequency of 5GHz ( Kundu et al.lll97^ 
and near the sol ar polar at 17GHz (Selhorst et al. 20inll . The result of higher frequency range (> 30GHz) 


(|2011r i. Please note that the limb brightening of the numerical results is 
higher than in observations. This may be attributed to the spicules, which will decrease the limb brightening 
(lElzneiill976r) . Another reason is perhaps no consideration of fluctuation, scattering effect, and magnetic 
field. We need more work and m ore observations to understand this problem. The new construction of 
the Chinese solar radio heliograph (jYan et al.l 120091 1 will provide more observations on the 2D image of the 
Sun at a wide frequency range of AQOMHz — IbGHz. In a word, this work gives an appropriate numerical 
method on the study of the quiet-Sun radio emission. 
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A. Appendix 


Here we discuss the the convergence of the integration in equation (j22[) . At the turning point poj 
/io = d/po) Mo is the refractive index at p = po- For p = po -I- d p, when dp — >■ 0, it can be proved that 
JL-^ -i 0. The refractive index y? = 1 — ( Benz 1993 1. where N = Nq dN cm.~^, Nq is the 
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electron density at p = po- So, = 1 — Let the constant Ci = 




. Thus, 


N^dp 


N^dp 


(Al) 


x/m' - y'l-Ci(iVo + diV)-6Vp2 

Presuming dA^ > 0 when dp > 0, we have ^ and Nq + dN < Nq. Thus Al can be transformed as 

N^dp N^dp 


< 


V^l-Ci(7Vo+dA^)-6Vp2 VI - CiNo - b^p-^ 

As 1 - CiNo - ^ = 0) and ^ = (p„+dp )2 = - (poPlpt+dpl^^' 

N'^dp N'^dp po{po + dp)N^ ^/d^ 


(A2) 


When dp —>■ 0, po + dp —>■ po, and A^ —>■ A^o- So 


dx/2po + dp 


(A3) 


A^2dp 




^ Po(po + dp)N'^y/dp ^ po(2po + dp)N'^^/dp ^ Po(A^o)V2podp 


dx/2po + dp 


dy/lpQ + dp 


0 


(A4) 


That is, the integration of equation (1221) near the turning point (singularity point) is convergent so long as 
the electron density A^ increased limitedly and monotonously along the height. 
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